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Abstract 

A precise determination of the charm-strange baryon lifetime is presented. The 
data were accumulated by the Fermilab high-energy photoproduction experiment 
FOCUS. The measurement is made with 300 H+ — ► H~7r + 7r + decays, 130 — > 
T, + K~n + decays, 45 — » pK~ir + decays and 58 — > A°A' _ 7r + 7r + decays. The 
H+ lifetime is measured to be 0.439±0.022±0.009 ps. 



1 Introduction 

The lifetime hierarchy of the weakly decaying charm mesons is well estab- 
lished [1]. However, the pattern of the predicted lifetimes for the weakly de- 
caying charm baryons agrees only qualitatively with experimental results [2,3]. 
The A+ lifetime is known to an accuracy of ~5% [4,5], but the others, 
S°, have uncertainties on the order of 20%. Interestingly, baryon sec- 
tor lifetime measurements provide information on quark interference and W- 
exchange. The essential difference from the mesons is that W-exchange among 
the valence quarks of the baryon is neither color nor helicity suppressed. The 
measured 5+ lifetime is larger than theory predicts, but there is a large ex- 
perimental uncertainty [6,7]. A more precise measurement could be conclusive 
in testing predictions in this sector. 

In the FOCUS spectrometer^ high energy photons with (E) « 180 GeV 
interact in a segmented BeO target to produce charmed particles. Charged 

1 FOCUS spectrometer is an upgraded version of the Fermilab E687 spectrome- 
ter [8]. 
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particles are tracked in the target region by two silicon vertex detectors which 
provide excellent vertex separation between the production and decay vertices. 
The average proper time resolution is ~ 50 fs for the modes used in this 
analysis. Downstream tracking and momentum measurement is performed by 
a system of five multiwire proportional chambers (MWPC) and two magnets 
of opposite polarity. The upstream magnet (Ml) is positioned downstream of 
the silicon detectors and in front of the MWPC system. The second magnet 
(M2) lies between the third and fourth MWPC stations. Charged particle 
identification is provided by three multicell threshold Cerenkov counters and 
two muon systems. One hadronic and two electromagnetic calorimeters are 
used to measure particle energy. 



2 Reconstruction of hyperons H and S + 

A detailed description of the S~ and S + reconstruction can be found elsewhere 
[9]. The S~ decays into A°7r~, with the A being reconstructed, if possible, 
through the decay A — > £>7r~.[^] The H~'s are reconstructed according to 
where and how the decay occurs. We consider four categories: Typel where 
the S~ decays promptly without leaving a track in the silicon detectors; Type2 
where the S~ decays after passing through the silicon and the A daughter 
is fully reconstructed; MV decays, which are topologically identical to Type2 
decays, but are reconstructed from the intersection of the three MWPC tracks 
from the H~ decay with the S~ silicon track; Kink decays where the H~ passes 
through the silicon and the A daughter is not reconstructed. If a S~ candidate 
is reconstructed both as Type2 and MV, then we choose the Type2 to avoid 
duplication. In Kink type decays, due to the failure to reconstruct the A , 
there is normally a two-fold ambiguity in the determination of the momentum. 
However, by requiring that the decay occurs in Ml we remove this ambiguity. 
Our data is composed of 9% Typel, 66% Type2, 6% MV and 19% Kink. We 
have more than one million reconstructed S~'s in all four types. 

The S + baryon can decay into pre or mr + . Each channel is studied separately. 
The reconstruction is similar to the Kink category described above. For the 
neutron in the S + — > nii + decay we require 0.3 < E/p < 2.0, where E is 
the energy deposited in the calorimeters and p is the reconstructed momen- 
tum. For the S + — ► pir decays, the protons must be positively identified in 
the Cerenkov system. As in the case of Kink type S~'s, we have a two- fold 
ambiguity in the determination of the S + momentum. This ambiguity can be 
partially removed for the E + — > niT + decay by using the location of energy 
deposition in the calorimeters. 

2 Throughout this paper, the charge conjugate state is implied whenever the decay 
mode of a particular charge is stated. 
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3 Reconstruction of 5+ candidates 



The S+ candidate is reconstructed using a candidate driven vertexing algo- 
rithm [8]. A candidate vertex is formed using the silicon track information 
of the decay daughters when available. This (secondary) vertex is required to 
have a confidence level (CL) above a value optimized for each topology re- 
constructed. We construct a seed track using the momentum vector of the H+ 
candidate and intersect it with at least two other tracks to form a primary ver- 
tex. This primary vertex is required to have a confidence level greater than 1% 
and to be in the target material to within 3 units of the error in the primary 
vertex position (TGM<3). We tighten this requirement to TGM<0 for the 
cases where the H~ or S + are partially reconstructed. Other cuts used to opti- 
mize the signal are the confidence level that any other track originates from the 
secondary vertex (IS02), the error in the proper time (a t ), and the significance 
of separation of the primary and secondary vertices (L/ctl). Cerenkov particle 
identification is accomplished by constructing x 2 like variables for the differ- 
ent particle hypotheses [10]. Briefly, we compute likelihoods for the various 
stable particle hypotheses e, ir, K and p. The pion consistency of a track is 
defined by a requirement on AW W = W min — W T , where W w is the negative log 
likelihood of the pion hypothesis and W m i n is the minimum negative log like- 
lihood of the other three hypotheses. Similarly, we define AWk = W n — Wk 
and AW P = — W p for use in identifying kaons and protons. We require 
AW W >-6 and AW K >2 for pion and kaon identification in all modes. Addi- 
tional particle identification is mode dependent and is described below. The 
resulting invariant mass distributions are shown in Figure 1; the combined 
plot is shown in Figure 4(a). 
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Fig. 1. E 7T + 7r + (4 reconstruction methods), Y, + K ir + (2 decay modes), pK tt^ 
and A i^~7r + 7r + invariant mass distributions with cuts as described in the text. 
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3.1 2+ — > S 7T + 7r + selection criteria 

The secondary vertex is required to have a CL greater than 2% except for 
Typel S~ where we require it to be greater than 10%. For Kink type decays 
we require that a t be less than 0.06 ps and that IS02 be less than 0.01%. The 
additional cuts are necessary since we cannot make an invariant mass cut on 
2~ candidates of this type. We demand that L/a L be greater than 4. 

3.2 S+ — > T l + K~n + selection criteria 

If there is an ambiguity in the S + momentum determination we resolve it as 
follows. First, we select the higher momentum solution when the difference 
in the masses calculated using each solution is less than 30 MeV/c 2 . When 
the difference is greater than 30 MeV/c 2 we reject a solution candidate if the 
other — > T l + K~n + mass is within 2.5a of the nominal mass. For a 
small fraction of cases this can lead to rejecting or accepting both solutions. 
Resolving the ambiguity in this way may create a bias and we have studied 
this as a possible source of systematic uncertainty. The CL of the secondary 
vertex must be greater than 2% and we require a t to be less than 0.07 ps. We 
ask that L/a L be greater than 7. 

3.3 S+ — > pK~n + selection criteria 

We identify protons by requiring that AW P > 10 and that the proton hy- 
pothesis is favored over the kaon hypothesis by two units of likelihood i.e. 
W K - W p > 2. We demand that W p - W K > for the kaon candidate. For 
pion candidates we require the momentum to be greater than 5 GeV/c. The 
secondary vertex must have a CL greater than 15% and IS02 less than 0.01%. 
To reduce combinatorial background, the candidate momentum is required 
to be less than 120 GeV/c. Long lifetime backgrounds from charmed mesons, 
where one particle is misidentified, are removed by making an invariant mass 
cut. In this way we eliminate contamination from the decays D + — > K~it + -n + , 
D + {Df) -> K-K+7V+ and D° -> K+K~. Finally, cuts of a t less than 0.08 ps 
and L/aL greater than 7 are applied. 

3.4 5+ — > A°A'~7r + 7r + selection criteria 

We select A candidates which decay downstream of the silicon vertex detector. 
To reduce contamination from K$ — > 7r + 7r~ decays, we require AW P > 8 for 
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the proton in the A decay. The secondary vertex must have a CL greater than 
15% and IS02 less than 0.01%. Cuts of a t less than 0.10 ps and L/a L greater 
than 4 are applied. 



4 Lifetime technique. 



We perform a binned maximum likelihood fit [11] to extract the lifetime of 
the S+ from the aforementioned decay channels. We fit the reduced proper 
time (£') distribution, defined as t' — (L — Na^)/ fl'-fc where N is the vertex 
detachment cut, (5c is the particle velocity and 7 is the Lorentz boost factor 
to the H+ center of mass frame. The t' distribution for S+ is of the form 
e~* / T , where r is the lifetime of the S+. A fit is made to the t' distribution 
for events which lie within ±2cx of the mass, where a is ~8 MeV/c 2 for 
the — > A°K~7r + 7T + channel and ~10 MeV/c 2 otherwise. The background 
is assumed to have the same lifetime behavior in the signal region and in the 
sidebands which are A-Yla away from the peak. Taking S as the number of 
signal events in the mass region and B as the total number of background 
events in the same region, the expected number of events n\ in the i th reduced 
proper time bin centered at t' is given by: 

rii = b — h d (1) 



where bi describes the background reduced proper time evolution as estimated 
from sidebands and /(^) is a correction function, which takes into account the 
effects of spectrometer acceptance, analysis cut efficiencies, and absorption of 
the particles as a function of the reduced proper time. In Figure 2 we plot 
the correction function for each decay mode. The large corrections in some 
modes at low reduced proper time are due to the suppression of short lived 
decays by our selection cuts. The likelihood is constructed from the product 
of the Poisson probability of observing events when rii are expected with 
the Poisson probability of observing X) bi background events when AB are 
expected. The factor of 4 accounts for the fact that the sideband region is four 
times wider than the signal region. The likelihood takes the form: 

^ ( n<r!i) X( (^^!) (2) 
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Fig. 2. Lifetime correction functions for the different decay modes and topologies 
described earlier. 

The combined likelihood function is given by the product of the eight likeli- 
hoods as shown in Eq. (3). 



fZL 7T' TV * 

'-'Typel 



X L "Type2 



X ''-E+(n7r+) X ^E+(p7r°) X *->pK-Tr+ x ^AOK-tt+tH 



(3) 



There are nine parameters in the fit, one parameter for the lifetime r and 
eight parameters for the backgrounds one for each decay type. Our result for 
the S+ lifetime with statistical errors is 0.435 ± 0.022 ps. 



5 Systematic Studies 



We compute the lifetime using several different L/cl cuts. The results are 
shown in Figure 3(a). Systematic effects were studied by computing the life- 
times of data samples split by individual H~ topologies and modes. The results 
are shown in Figure 3(b). All variations are consistent within statistical un- 
certainties and do not contribute to a systematic uncertainty. We have inves- 
tigated systematic effects due to the t' resolution by examining the variance 
in the fitted lifetime for different if bin size, by reducing the t' range for the 
fit from 3 to 2 ps and by excluding the lowest t' bin from the fit. Uncertain- 
ties in the measurement of particle momenta can lead to a systematic shift 
in the reduced proper time. Our studies show this shift to be 4 ± 2 fs. The 
final quoted value is adjusted by this amount and a systematic uncertainty of 
2 fs is included in the final systematic uncertainty. The treatment of the two 
solution ambiguity in the T l + K~n + mode creates a small bias in the measured 
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Fig. 3. a) Lifetime stability versus L/cr^; in parenthesis is the L/ct,l cut used for 
T< + K~it + and p.KT~7r + . b) lifetime measurements for systematic studies. The solid 
line represent the central value with the dotted lines showing the extent of the 
statistical error at L/cjl > 4(7). 

Table 1 

Contributions to the systematic uncertainty. 



Contribution 


Uncertainty (fs) 


Momentum 


±2 


two solution bias 


±1 


Split sample 





t' Resolution 


±4 


Background 


±2 


Fit variant 


±8 


Total 


±9 



lifetime due to an overestimation of the background in the signal region. This 
effect is less than 1 fs on the total lifetime. Systematic effects due to the back- 
ground were investigated by varying the width of the sideband regions, and by 
altering the background level by imposing a minimum separation between the 
primary and secondary vertices of 1.5 mm. The variance from these tests is 
added to the systematic uncertainty. We tested different fit conditions; exclud- 
ing MV type decays from the lifetime fit, taking the weighted average of the 
split samples, and using a combined fit to the mass shape and reduced proper 
time of — > £ + A'~7r + as an alternative method of treating the two solution 
ambiguity. The systematic uncertainty due to the variation in fit conditions is 
taken to be the sample variance since we consider all of the measurements to 



8 




Fig. 4. a) The invariant mass for the combined sample, b) The combined lifetime 
fit of the H+ modes with a background subtracted, Monte Carlo corrected, reduced 
proper time distribution. 

be equally valid. The systematic contribution is found to be 8 fs. 

The components of the systematic error are presented in Table 1. Adding these 
contribution in quadrature gives a total systematic uncertainty of 9 fs. Figure 
4(b) shows the background subtracted, Monte Carlo corrected, reduced proper 
time distribution for the S+ signal. 



6 Summary 



We have measured the using four decay modes which occur in eight dif- 
ferent topologies with a combined sample of 532.4 ± 30.4 events. We find the 
lifetime to be 0.439 ±0.022 ±0.009 ps where the first error is statistical and the 
second is systematic. Our result is compared with other experimental values 
in Table 2. Our measurement uncertainty is a factor of two better than those 
of the current world average [1]. As discussed in the introduction, theoretical 
models predict the lifetime hierarchy of the charmed baryons. A number of 
authors [2,15-17] predict that t(S+) > t(A+) where the inequality represents 
a factor of about 1.3. Using the A+ lifetime average of PDG, CLEO, and SE- 
LEX [1,4,5] (0.1916 ± 0.0054 ps) and the H+ lifetime reported in this paper, 
a ratio t(S+)/t(A+) = 2.29 ± 0.14, is obtained. Our well measured ratio is 
significantly different from predictions of order 1.3. 
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Table 2 

lifetime measurements 



Experiment 


Lifetime (ps) 


Events 


Year 


WA62 [12] 


n A o+0.21+0.20 
U - 48 -0.15-0.15 


53 


1985 


E400 [13] 


0.40l§;J| ±0.10 


102 


1987 


ACCMOR (NA32)[14] 




6 


1989 


E687 [7] 


0.411^ ±0.02 


30 


1993 


E687 [6] 


0.34±g:g| ± 0.02 


56 


1998 


PDG (average) [1] 


n qq+0.06 




2000 


This Measurement 


0.439 ± 0.022 ± 0.009 


532 


2001 
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